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I. INTRODUCTION 



Quantum dots(QDs) are considered to be promising as the elementary basis for the new 
generation of semiconductor devices [JJ [2] . The unique opportunity to perform the energy 
level control and flexible manipulation in QDs is due to the full quantization of charge carrier 
energy spectra in these systems. This allows design and manufacturing of artificial structures 
with prescribed quantum physical characteristics [3] . That is why the scope of QDs potential 
applications is very wide, from heterostructure lasers to nanomedicine and nanobiology. An 
impressive example of such application is represented by QD lasers possessing low threshold 
current and high efficiency [3j- 

The peculiarities of physical processes in QDs are caused by both their composition and 
geometry. Electronic, kinetic, optical and other properties of QDs have been investigated 
experimentally and theoretically in many papers [H-T13]. Particularly, the optical absorption 
characteristics of QDs have been shown to be strongly correlated with their geometry, on 
one hand, and with their physical-chemical properties, on the other hand. In one of the 
first publications on optical transitions in QD [T3] the interband absorption of light was 
considered in the ensemble of weakly interacting spherical QDs implanted in a dielectric 
matrix. The dispersion of QD sizes was characterized in the framework of Lifshitz-Slezov 
theory [15] . It was shown that in the absence of size dispersion, due to the full quantization of 
charge carriers energy spectra in QD, the absorption coefficient behaves like a delta function, 
and the absorption threshold frequencies depend on the peculiarities of electron and hole 
energy spectra. When the QD size dispersion is taken into account, the averaging procedure 
yields the absorption profile having finite width and height. 

Recently several reports concerning the experimental implementation of narrow-band 
InSb QDs have appeared [T6l [T7] , in which the dispersion law for electrons and light holes 
is non-parabolic and described according to the double-band mirror Kane model [TBI HH]- 
For non-interacting band of heavy holes the dispersion law is considered as quadratic. The 
investigation of optical absorption peculiarities in InSb QDs with the transitions from light 
and heavy hole bands to the conduction band taken into account is an interesting problem. 
Interband transitions in an ensemble of cylindrical or spherical InSb QDs were considered 
theoretically in the dipole approximation with and without magnetic field, including exciton 
effects, by means of the perturbation theory and the adiabatic methods j20H22]- In our ear- 
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lier work we elaborated the calculation schemes, symbolic-numerical algorithms (SNAs) and 
programs, based on the generalized Kantorovich method (KM) for numerical solving with 
required accuracy the boundary-value problems (BVPs) of discrete and continuous spectra 
describing the axial-symmetric models of quantum wells(QWs), quantum wires(QWrs) and 
quantum dots(QDs) in external fields within the framework of the effective mass approxi- 
mation [23H35] . Meanwhile, for the analysis and estimations of the appropriate range of 
material parameters, spectral and optical characteristic of quantum dots at the first stage of 
investigation, conventionally, approximate eigenvalues and eigenfunctions evaluated in the 
analytical form were applied [BHE1 EEH 122] • However, it is a real challenge to specify the range 
of applicability of such approximations in the problems, depending on a few parameters [2] , 
e.g., for impurity states of quantum wires in a homogeneous magnetic field [25J. 

With this aim in the present paper we report the formulation and MAPLE- environment 
implementation of algebraic schemes of the perturbation theory (PT) of the Lennard- Jones 
(LJ) and Rayleigh-Schrodinger (RS) [30], permissive in the nondiagonal and diagonal adia- 
batic approximations, respectively, to evaluate in numerical and in analytic forms the eigen- 
values and eigenfunctions of models of spheroidal QDs in homogeneous magnetic and electric 
fields. To construct the required perturbation schemes, we choose such models of spheroidal 
QDs, in which the basis functions depending upon fast variables can be expressed in the 
analytic form. The region of the model parameters, for which the PT asymptotic series are 
applied, is estimated using the results of numerical calculations carried out with required 
accuracy. The efficiency of the schemes is demonstrated by the analysis of spectral charac- 
teristics of oblate and prolate spheroidal QDs and also spherical QDs with corresponding 
shape of confinement well with walls of infinite height under the influence of homogeneous 
electric fields (HEFs). We apply the developed approach to the analysis of spectral charac- 
teristics of oblate and prolate spheroidal QDs with parabolic and non-parabolic dispersion 
laws under the influence of HEFs, i.e., the quantum-confined Stark effect. 

The paper is organized as follows. In Section 2 the calculation scheme for solving elliptic 
BVP describing spheroidal QDs in homogeneous electric fields using the Kantorovich method 
is presented. Section 3 is devoted to the description of the PT schemes by slow variables 
in nondiagonal adiabatic approximation and the comparison of the results with those of 
numerical calculation with given accuracy. In section 4 the explicit PT scheme for evaluation 
of the basis functions of the fast variable for oblate spheroidal QDs in a homogeneous electric 
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field is derived. Section 5 is devoted to the description of PT schemes by slow variables in the 
diagonal adiabatic approximation for spheroidal QDs in electric fields. The results evaluated 
here in the analytic form are compared with numerical ones to establish the range of their 
applicability. In Section 6 the absorption coefficient for an ensemble of spheroidal QDs with 
random dimensions of minor semiaxis and with parabolic and non-parabolic dispersion laws 
for holes and electrons under the influence of HEFs is found using the calculated eigenvalues 
and eigenfunctions. In conclusion we summarize the results and discuss further applications. 



II. STATEMENT OF THE PROBLEM 



Let us consider an impurity localized in the center of a quantum dot and take the electron- 
hole interaction into account. Then in the effective mass approximation of the k • p-theory 
the Schrodinger equation for the slow- varying envelope wave function ^(f e , f/J of an electron 
(e) and a hole (h) in a uniform magnetic field H with the vector-potential A = x f and 
electric field F in oblate and prolate QDs reads as [8]: 



{H(f e ,f h )-E}v(f e ,f h ) = 0, (1) 
H(f e ,f h ) = W-Lfe- ^ A ) 2 - qi (F • r~) + U mf {fi) - hft) + 

i=e,h v 



QeQh 



Here fj is the radius- vector, |fj| = \/ x 2 + + z 2 , p, = -«W r - is the momentum, E is the 
energy of the particles, q e = — e, q^ = +e, and q c are the Coulomb charges of the electron, 
the hole, and the impurity center, k is the dc permittivity, \ii = (3 e (h)mo is the effective mass 
of electron or hole, mo is the mass of electron. For the model under consideration, U(r) is 
the potential of a spherical or axially-symmetric well 

U(fi) = {0, S(fi) < 0; L> , S(f t ) > 0}, (2) 

bounded by the surface S(fi) = with walls of infinite height (infinite potential barrier 
model, IPBM) or finite height 1 Uq < oo (finite potential barrier model, FPBM). In Eq. 
(J2j) S(fi) depends on the parameters a, c, which are semiaxes of a spheroidal QD, 

5(f 4 )^(x 2 + y 2 )/a 2 + .t/c 2 -l. (3) 
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Below we restrict ourselves to IPBMs of spheroidal quantum dots with possible influence 
of the uniform electric field F = (0,0, F), the magnetic field being switched off, H = 0, 
and the Coulomb interaction of the electron and the hole with the impurity center being 
absent, q c = 0. In this case the wave function ^(f e ,fh) = ^ e (f e )^ h (fh) is factorized. So, 
we arrive at the 3D BVPs for unknowns \£ /6 (r e ) and E e or ty h (fh) and E h . The eigenvalues 
and eigenfunctions needed to evaluate the absorption coefficients (ACs) were calculated with 
prescribed accuracy by means of the program packages ODPEVP and KANTBP [28-30J. 
The models with nonzero values of these parameters were announced in [HI [2S] . Throughout 
the paper we make use of the reduced atomic units [21 13]: a* B = nh 2 / p p e 2 is the reduced 
Bohr radius, Er = Ry* = h 2 /(2p p a* B 2 ) is the reduced Rydberg unit of energy, and the 
following dimensionless quantities are introduced: *(f) = a B ~ 3/2 ^(r), 2H = H/Ry% S = 
2E = E/Ry\ 2U(r) = U(r)/Ry*, r = r/a* B , a = a/a* B , c = c/a* B , 2 lF = F/Fj, F* = 
Ry*/(ea* B ) = e/(2 K (a B ) 2 ). 



A. The BVP for SQDs in the effective mass approximation 

In cylindrical coordinates z,p,cp the solution of Eq. Q, periodical with respect to the 
azimuthal angle (p, is sought in the form of a product ^(p, z, </?) = \l/ m (p, 2;)exp (imp) / y/2n , 
where m — 0, ±1, ±2, ... is the magnetic quantum number. The 3D BVP for SQDs at fixed 
values of m is reduced to 2D BVP with respect to fast Xf and slow x s variables: oblate 
Xf = z (minor axis), x s = p (major axis) and prolate Xf — p (minor axis), x s = z (major 
axis) [27]: 

(Hf(x f ;x a ) + H s (x s ) + V fs (x fl x s ) -£ t m ) *rO&/,z«) = °- ( 4 ) 

Here H s (x s ) is the operator of slow subsystem 

Id d 

H s (x s ) = — r— g 2 s(x a )n V s( x s), (5) 

gu(x s )dx s dx s 

and Hf(xf] x s ) is the operator of fast subsystem 

H f (xf, x a ) = + V f (x f] x s ). (6) 

For OSQD g ls (x s ) = g 2s (x s ) = 1, gif(x s ) = 92f(x s ) = p, V f (x f ;x s ) = 0, V s (x s ) = m 2 /p 2 , 
Vf s (x f ,x s ) = 2j F z, while for PSQD g ls (x s ) = g 2s (x s ) = p, gi f (x s ) = g 2 f{x s ) = 1, 



Vf(xf,x s ) = m 2 /p 2 , V a (x s ) = 2^ F z, Vf s (xf,x s ) = 0. From (|2j) the boundary conditions 
for the eigenfunctions ty™(xf, x s ) of SQDs, corresponding to a well with walls of infinite 
height, have the form 



lim p- 

p->0 1 



dp 



{p,z} 



-5 0m + ^(p,z)(l-5 l 



Om, 



0. 



< 1 



{p,z} 



p 2 z 2 

— H 

a 2 c 2 



ao 2 
= 1 



The eigenfunctions *i?™(xf,x s ) corresponding to the eigenvalues £™ = E™ < £™,... are 
subject to the normalization and orthogonality conditions 



/ pdpdz*?(p,z)*?(p,z) = 6 t r. 
Jn 2 



Note, that at 7^ = the solutions are separated by the z-parity er = ±1 into two invariant 
subspaces ^J™ 7 corresponding to the eigenvalues £™ G = £™° < £™ a , ... , while at 7^ 7^ the 
z-parity is broken. 



B. Kantorovich or adiabatic reduction of the BVP 



The solution ty t (xf,x s ) = ^^{xf, x s ) of the above problem at fixed m is sought in the 
form of Kantorovich expansion 

J max 

%(x f ,x a ) = ^B j {xf,x s )x jt {xs)- (7) 
i=i 

The set of appropriate trial functions is chosen as the set of eigenfunctions Bj(xf,x s ) cor- 
responding to the eigenvalues Ej(x s ) of the Hamiltonian Hf(xf;x s ), Eq. (JHl, depending 
parametrically on x s G VL{x s ): 

H f (x f ;x s )B j (xf,x s ) = E j (x s )B j (xf 1 x s ). 

The eigenfunctions Bjxf, x s ) corresponding to the eigenvalues EAx s ) = E^Xg) < E 2 (x s ), ... 
are subject to the normalization and orthogonality conditions with the weighting function 
gif{xf) in the same interval Xf G Q Xf (x s ): 

/ B i (xf,x s )B j (xf,x s )g lf (x f )dxf = 5 i:j . (8) 

Jxf in (x 3 ) 
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The BVP for a set of ODEs of the slow subsystem with respect to the unknown vector func- 
tions Xti x s) — (Xi;t( x s), ■■■■> Xj m ^;t( x s)) T corresponding to the unknown eigenvalues 2E t = £ t , 

(d + E(x s ) + W(x s ) - IS^j Xt ( x s) = 0, 



wrf \ t t/ \ I 92s{x s ) 1 dg 2s (x s )Q(x s ) g 2s (x s ) d 

W(x s ) = U{x s ) + — —H(x s ) + — -— + — r^Q( Xs )^~ 

9is{x s ) gis{x s ) dx s 9is{ x s) dx s 

satisfy the orthogonality and normalization conditions 

/ (Xt( x s)) T Xt'( x s)9is(x s )dx s = 8tf. (10) 



pxf ax (x s ) 

Uij(x s ) = Uji(x s ) = / B i (x f ;x s )V fs (x f ,x s )B j (x f ;x s )g lf (x f )dx f , 

Jxf in (x 3 ) 



Here the effective potentials Hij(x s ) and Qij(x s ) are defined by the formula 

H.Ax,) = Hji ( X .) = r" 1 "' dB ' { ^ X - )dB ^''- X - ) 9lf (x 1 )±c f , (11) 
Qij{x s ) = -Qji(x s ) = - Bi(x f ;x s ) — 3 f ' s gif(x f )dxf. 

Jxf in (x 3 ) ax s 

Here the basis functions of the fast subsystem and the matrix elements are calculated ana- 
lytically. For oblate spheroidal QDs (xf = z, x s = p) 



B i (x f ; x s ) = B° (x f ; x s )= sin 



ixn 



c-sja 2 — x 2 s \ 2 \ ca/1 — x 2 s /a 2 



a 2 „ n 2 i 2 



E i (x 8 ) = E?(x,) = E ii0 E v ^ = —, U ii (x a )=0, (12) 

U(x)-U (x) ^^ U (x) _ ^Fcij(-l + (-iy+>) 

a K 1, J ) " 

Hu{x s ) = Hu;o(x s ) ^ 2 — ^ 2 , Hu;o(x s ) = ^2 a 2 ' 

h(x)-h (x) q2 ^ h (x) _ m 2 +m+(-iy +j ) 

n l j\x s ) — si lj;0 {x s ) ^ 2 _ 2 , n lj]0 yx s ) - a 2 (i 2 - j 2 ) 2 ' 

n< \-n ( \ aXs n i \ - Mi±lzll!!!) ■ +• 

^ij\ x s) — ^ij;0{X s ) ^ 2 _ rt-2 1 tyij;0\ x s) — ( '2 _ '2\ ' J I % - 

a\% — J ) 
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For prolate spheroidal QDs (xf = p,x s = z) (at m = for nondiagonal potentials i ^ j) 



- 2 («/' ) 2 

= Ei-oj-^——, E i]0 = ™l , (13) 



(cr — x^J a z 



£/n(rr s ) = 0, C/y(x s ) = 0, 

, s _ , > C 2 X 2 I 1 + (^m|) 2 ) 

±lii\X s ) Ilii-QyX s ) . 2 2^2 ' ii;Q a y q 2 ' 

rr/ N_rr / X C ^ rr / x _ 2 / = jj Jl(Jfc) J^x) , 

H„(.r,) - /iyjQ^.J ^ 2 _ ^ 2 , ily-o^aj - ^ I -Vo / Jl(4) Ji(J^ X 



<r — 



J n I _ " y — * x — ■ J n I _ > — ■ x — ■ ^ 

°io Ji(J3) Ji(Jg) ^0 ./i(Jg) 

Qij{Xs) = Qij;0{X s )-^ -3, Qij-o{X s ) = -^0 / J ( Ji\ j / jj\ X ® X i J T % i 

c —x s C J Jl{J ) Ji{J J ) 

where a npp+ i : \ m \ = J^f +1 are positive zeros of the Bessel function of the first kind [37]. 

For the interesting lower part of the spectrum S t : S\ < S 2 < the number j max of 
the equations solved should be at least not less than the number of the energy levels of 
the problem ([9j at a — c — tq. To ensure the prescribed accuracy of calculation of the 
lower part of the spectrum discussed below with eight significant digits we used j max = 16 
basis functions in the expansion (8) and the discrete approximation of the desired solution 
by Lagrange finite elements of the fourth order with respect to the grid pitch £1^ (x s ) = 
[x S ;mm, %s;k = %s;k-i + h s ; x s;max ]. The details of the corresponding computational scheme are 
given in [24|. 

The convergence of eigenenergies St vs number j max of basis functions for oblate and 
prolate spheroidal QDs, and for spherical QD is shown on Tables |T] and [II] at = and 
m = 0. The considered QDs having the size comparable with De Broglie wavelength of 
composed particles with small effective masses are referred as quantum-size systems. In the 
spheroidal QDs having different length of minor and major axes the quantization procedure 
leads to different transversal and longitudinal spectra. Moreover, for PSQD (c = 2.5, a = 
0.5) the confinement in two variables (xy) with the minor semiaxis a = 0.5 leads to greater 
eigenvalues, than the confinement in one variable (z) with the size-for-size minor semiaxis a = 
0.5 for PSQD (c = 2.5, a = 0.5). Tables [I] and [lT| show that the expansions in basis functions 



(12) and (13) in cylindrical coordinates have better rate of convergence in the adiabatic 
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TABLE I: The convergence of eigenenergy St vs number j max of basis functions at 7^ = . 
Fast and slow variables Xf = z and x s = p (oblate SQD and spherical QD), number of nodes 
i = (n-zo = n a — 1, n po ), * notes diagonal approximation at j = 2 





a = 2.5, c = 0.5 


a = 2.5, c = 2.5 


I Tljyr, . Tb/jn ) 


(0,0) 


(0,1) 


(2,0) 


(0,0) 


(0,1) 


(2,0) 


c 


12.737 41 


19.936 21 


96.696 83* 


1.468 496 


5.445 665* 


5.589 461 


1 


12.765 48 


20.046 02 


96.753 17* 


1.590 238 


5.766 612* 


6.004 794 


2 


12.764 90 


20.041 33 


96.754 27 


1.580 243 


5.340 214 


6.329 334 


4 


12.764 82 


20.040 74 


96.752 15 


1.579 273 


5.316 872 


6.317 204 


16 


12.764 81 


20.040 65 


96.752 01 


1.579 140 


5.314 832 


6.316 562 


Exact 








1.579 136 


5.314 793 


6.316 546 



TABLE II: The convergence of eigenenergy £ t vs number j max of basis functions at 77? = . 
Fast and slow variables Xf = p and x s = z (prolate SQD and spherical QD), number of nodes 
i = ( n ppi n zp ), * notes diagonal approximation at j = 2 



Jmax 


c = 2.5, a = 0.5 


c = 2.5, a = 2.5 




(0,0) 


(0,2) 


(1,0) 


(0,0) 


(0,2) 


(1,0) 


c 


25.184 73 


34.428 85 


126.424 5* 


1.493 612 


5.131 784 


5.898 668* 


1 


25.201 74 


34.530 30 


126.456 5* 


1.584 433 


5.680 831 


6.071 435* 


2 


25.201 29 


34.525 78 


126.457 3 


1.579 860 


5.331 101 


6.324 717 


4 


25.201 21 


34.525 12 


126.456 1 


1.579 239 


5.316 732 


6.317 058 


16 


25.201 20 


34.525 02 


126.456 1 


1.579 138 


5.314 828 


6.316 554 


Exact 








1.579 136 


5.314 793 


6.316 546 



limit of strongly oblate and prolate QDS, than for the benchmark spherical QDs with the 
known spectrum, which is not surprising. For lower states the crude adiabatic approximation 
(without Hjj(x s )) (CAA) provides a lower estimate, while the adiabatic approximation (AA) 
(with Hjj(x s )) (1) gives an upper estimate, such that at the ratio of minor to major semiaxis 
equal to 1/5 the bracket is approximated with the accuracy of ~ 0.1%. 

Below we present the analysis of the spectrum under the variation of parameters, which 
opens the questions about the additional symmetry of the problem, connected with the 
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existence of exact and approximate integrals of motion [271 1381 




FIG. 1: Dependence of eigenenergies £ (in units of E e ) of lower part of spectrum of electronic 
states of QDs at m = on electric field strength 7^ (in units of Fq): for spherical quantum dot 
(SQD) with radius a = c = 2.5, oblate and prolate spheroidal quantum dots (OSQD and PSQD) at 
different minor semiaxis (for OSQD c = 0.5, 1, 1.5, 2, a = 2.5, for PSQD c = 2.5, a = 0.5, 1, 1.5, 2). 

In Fig. [I] we show the eigenenergies of the lower part of the spectrum £ t , t = 1, ...,40 at 
m = for OSQD (c = 0.5, 1, 1.5, 2, a = 2.5), SQD (c = 2.5, a = 2.5) and PSQD (c = 2.5, a = 
0.5, 1, 1.5, 2) as functions of the dimensionless strength 7^ of the electric field. In spite of 
the fact that at 7^ = the eigenfunctions of SQD, OSQD and PSQD have definite z-parity, 
and, therefore, exhibit additional integrals of motion and separation of variables in spherical 
and spheroidal coordinates systems, the spectrum of eigenvalues at fixed m is simple, i.e., 
nondegenerate, similar to the case 7^ 7^ 0, when the eigenfunctions have no definite z-parity. 
At 7^ = a one-to-one correspondence rule n pp + 1 = n p = i = n = n r + 1, i = 1, 2, ... 
and n zp = I — \m\ holds between the quantum numbers (n, l,m,a = (— l)' m 'cr) of SQD with 
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the radius r Q = a = c, the spheroidal quantum numbers {n^ = n r ,n v = I — \m\,m, a} 
of PSQD with the major c and the minor a semiaxes, and the adiabatic set of quantum 



numbers \n r 



n 



+ 1, n zp , m, a] under the continuous variation of the parameter ( ac = a/c. 



At 7jt = there is a one-to-one correspondence rule n Q = n zo + 1 



2n 



n = 1,2,3, 



and n po = (I — \m\ — (1 — a)/2)/2, between the sets of spherical quantum 



numbers (n,l,m,a = (— l)' m 'cr) of SQD with the radius r$ 



c and spheroidal ones 



I — \m\,m, cr} of OSQD with the major a and the minor c semiaxes, and the 



adiabatic set of cylindrical quantum numbers [n = n zo + 1, n po , m, a] under the continuous 
variation of the parameter ( ca = c/a. 
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FIG. 2: Eigenfunctions of sixth order of PT of 2D BVP for oblate SQD a = 5, c = 0.5, [t = n zo = 0, 
n = n Q = 1,2,3,4, m = 0] in electric field 717 = —10 (weak asymmetry by z - axis i.e. by 
minor ellipsoid axis). Eigenfunctions of sixth order of PT of 2D BVP for prolate SQD c = 2.5, 
a = 0.5, [t = n pp = 0, n = n p = 0, 1, 2, 3, 4, m = 0] in electric field 717 = —1 (asymmetry by z - axis 
i.e. by major ellipsoid axis). 

One can see that when the parameter 717 increases the eigenvalues St decrease faster for 
SQD, slower for PSQD and even more slower for OSQD, because the influence of the electric 
field for OSQD at c = 0.5 is essentially weaker than for PSQD at c = 2.5. With increasing 

a series of exact crossings of eigenenergies with different values of quantum numbers 
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for PSQD and OSQD occur at > 20 and a series of avoided crossings for SQD occur 
at 7^ > 10. With further growth of the parameter they first increase and then begin to 
decrease. Indeed, with the growth of 7^ the eigenfunctions with smaller number of nodes 
in the longitudinal variable z are localized (see Figj2| in the vicinity of the equilibrium 
point, and the corresponding eigenenergies decrease. Increasing the number of nodes is 
accompanied with derealization of the wave functions, and the corresponding eigenenergies 
increase and then decrease again. For PSQD the density of states per unit energy for the 
eigenfunction with the same number of nodes n pp in the transverse variable p is greater (i.e., 
the separation between the adjacent energy levels is smaller) than the density of states for 
the function having the same number of nodes n zp in the longitudinal variable z. For this 
reason in Fig. [I] one can see three crossing series of curves with different number of p-nodes 
n pp = 0,1,2, the lower of them (e.g., with a = 0.5, n pp = 0, and n zp from to 12) are 
decreasing at all jp > 0, while the upper ones (e.g., with a = 0.5, n pp = 0, and n zp starting 
from 13) with the energies, exceeding that of the state (n pp = 1) without z-nodes (n zp = 0), 
increase from the beginning and then start to decrease. Thus, at small 7^ the energy levels 
for the groups of states with even n pp = 0,2... and odd n pp = 1,3... number of nodes are 
repulsing and crossing. 

For OSQD, on the contrary, the number of energy levels per unit energy for the eigen- 
functions having the same number n po of p-nodes is smaller (i.e., the separation between 
the adjacent levels is larger) than that for the eigenfunctions having the same number n zo 
of z-nodes. Therefore, in Fig. [I] one can see four crossing series of almost 'parallel' curves 
with different number n zo = 0, 1, 2, 3 of z-nodes. 

For OSQD and PSQD the crossings of the energy levels that occur with increasing 7^ are 
similar to the exact crossings of the energy levels with decreasing c semiaxis in OSQD and 
PSQD without electric field (7^ = 0), i.e., we observe the accidental degeneracy, which is 
known to be generally associated with the existence of an additional integral of motion [2"T] 
and with the separability of variables in oblate and prolate spheroidal coordinate systems. 
Thus, from our observations it follows that an additional approximate integral of motion 
should exist. 

For SQD eigenfunction with different numbers of p- and z-nodes, n p and n z , and with 
increasing jp the series of crossings become mixed. Note, that the eigenenergies of the states 
with the same z-parity at 7^ = are repulsed with increasing 7^ (e.g., [t = 9, n — 1, I — 5, 
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£ 9 ( 7F = 0) = 14.01] and [t = 10, n = 3, I = 0, £ w (-fF = 0) = 14.21]), but the states with 
different z-parity are attracted (e.g., [t — 7, n — 1, I — 4, £7(7^ = 0) = 10.71] and [t = 8, 
n = 2, / = 2, £s(lF = 0) = 13.24]). This fact should be also associated with the existence of 
approximate integrals of motion. Indeed, from Fig[T]one can see that for SQD at a = c = 2.5 
with increasing the series of exact crossings appear. 



III. THE PTLJ IN NONDIAGONAL ADIABATIC APPROXIMATION 



We expand the potentials Q and Q of the BVP @ and (|To[ in Taylor series in the 



vicinity of x s = 0: 

^max 

Ei(x s ) = Ei, + Y^ 



Ei-0 ok 



-x c 



k=l 



Uij(x s ) = U ij;0 + ^2 



max j j 

U ij;k 2k 
X c , 



-2k 



(14) 



k=l 



Eij (x s 



n-max t j ^rnax s \ 

El ffijjO 2k n / \ _ ^k-\ 



k=l 



k=l 



T 



2k-l°°s 



where Uij-k = ^j^rUij-o and the parameter r equals r = a for OSQD, and r = c for PSQD. 
Substitution of expansions Q into Eq. Q leads to the BVP for a set of ODEs of slow 
subsystem with respect to the unknown vector functions Xt( x s) = (Xi;t( x s), •••> Xj ma .^,t( x s)) T 
corresponding to the unknown eigenvalues 2E t = £ t : 



Xi;t{Xs) 



Jmax fcmax 

i^j fc=i 



U if,k 2k , j H ij;0 2k , / 9 . _ 1 \ Q»j;0 2fc-2 



«i;0 2fc 



-2A,-l X s 



(15) 



0. 



where E/y.fc is given by the expansion (14) and V s (x s ) = for OSDQ; Uij(x s ) = and 
V a (x s ) = ^fz for PSDQ. We choose the unperturbed operator to have the eigenvalues and 
basis functions of 2D and ID oscillators. For the OSQD (2D oscillator) with respect to the 



scaled slow variable x we have: x s = p = J x/ y/Ef), where Ef = (Ei>. + Hj/ i /. )/(4 
u;?//4, i.e., the adiabatic frequency, at given i' = n Q 



L(n)=V^-E^\ D <°> = -(^-f-£), EM = E^ = n + (\m\ + l)/2, 

^\ X \rn\/2 { _ x/2)L \m\ {x) ,00 

T> , | m , / $ g = S qq/ . (16) 
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Therefore, the action of the operators L(n) and x on the function <S> q °\x) = (x) is 
determined by the recurrence relations 



L{n)*$l(x) = (q-n)*W n (x) 



■ q,m 

x$l°i(x 



q,m \ 



g+|m| v ^$f_ ) lim (x) + 



+ (2g + \m\ + 1)$® (z) - v/ g +H + lv^Tl$S lim (x 



(17) 



For PSQD (ID oscillator) with respect to the scaled slow variable x x s = xj \/Ef, where 
Ef = (Ei'-o + Hiii'. )/c 2 = u 2 , i.e., the adiabatic frequency, at given %' = n p , we have 

d 2 



L(n) = n<®-E<-°\ D 



(o) 



dx 2 

exp(-x 2 /2)H q (x) 
^pK^2fl^fq\ 



+ x 2 , E (0) = E^ = 2n + 1, n = 0,l,...., 

/oo 
<$> q °\x)<$> q °\x)dx = 5 qq , (18) 
-oo 



Correspondingly action of operators L(n), x and 4- on function$g°' ) (x)is determined by 



recurrence relations 13 



L{n)^\x) = 2(q-n)$ q 0) (x) 



x <S>(0) (x) = ^ q %( x ) + ^S&g! (x) 



d 

c/x 9 



v/2 



(19) 



$l 0) fx) = - ^^^(x). 



V2 



V2 



The eigenfunctions (15) as functions of the new scaled variable x are sought in the form 
of expansion over the basis of the normalized functions $q (x), q — 0, 1, .... of the 2D or ID 



oscillators with unknown coefficients bj jS : 



XjA x ) = ^2 b j,i;t® g 0) (x), h, „. b 

q=0 



q<0;t — ^j,g>g max ;t 



0. 



(20) 



Below we demonstrate that such expansions are appropriate for getting approximate solu- 



tions in the lower part of the BVP spectrum (\9h and (10). Substitution of the expansion 
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(20) into (15) yields the set of equations 



<?max 



J max <?m ax 



9max 



(21) 



q=0 



jjki=l q=0 



q=0 



A a 2^ 1 2fcp (H-l)/2 x * + 

fe=l \ T 



2fc-2 , r,™2fc-l ^ 



(2fc - l)x^ + 2x 



where k = 2 and K(x s ) = for OSQD; k — 1 and V^(ar) = for PSQD. Applying the 



relations (17) or (19) to get first derivatives of the basis functions, we get the expressions 



for the action of operators A;,-: 



<?max 



q'=0 

and, hence, the algebraic eigenvalue problem with respect to the unknown E t and bj^ t 



(22) 



Qmax 



Jmax Qxnax 



&ii;q'qbi,q;t + Olij;q'qbj,q-t — K 2 ^tEj 1 2 &i,g;t. 

q=0 j^i=l q=0 



(23) 



In the matrix form it reads as 



AB t = K- 2 S t EJ 1/2 B t , Bj,B t = 8 tt >, 



where B t = (&i, 0; t, h, qmaxi ;t, h,o-,t, bj m ^, q ^-t) T is a vector with dimension of 

jmax(<7max + 1), and A is a positive defined symmetric matrix having the dimensions 
(jmax(g m ax + 1)) x (j max (g max + 1)) with the elements A iq ^ +m „ l)+q+h{qm ^ +1)(j _ l)+q , +1 = 

a ij;qq'- 

Note, that the approximation with nonzero elements on the diagonal of the matrix A = 



{ a iMq'qs q ^T=Q^i=io,o=io-i obtained by the action of the diagonal operator A^, Eq. (|2l|), on the 



basis function $$ 0) (x), Eq.(|22j), gives the diagonal adiabatic approximation (AA) of PTLJ 



solution (23), i.e., £ t ~ £i- n ,n = 0, 1, ... at each fixed i. Such adiabatic classification of the 



eigenenergies is used in Tables discussed below. 



The convergence of eigenenergies of Eq. (23) vs the order fc max of approximation of the 



effective potentials (14) for j max = 4 and g max = 60 is shown in Tables III and IV for OSDD, 



PSQD, and SQD at 7 F = and in Table M at j F = -10 for PSQD and SQD. Table IV 
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TABLE III: The convergence of eigenenergies £t of Eq. (23) vs order fc max of approximation of 
effective potentials from (14) for j max = 4 and g max = 60 basis functions at = 0. Fast and 
slow variables Xf = z and x s = p (oblate SQD and spherical QD), number of nodes i = (n zo = 

TL 1, Tlpo) . 



^max 


a = 2.5, c = 0.5 


a = 2.5, c = 2.5 


{jl Z oi IT'po) 


(0,0) 


(o,D 


(2,0) 


(0,0) 


(0,1) 


(2,0) 


8 

12 
20 


12.668 20 
12.749 67 
12.784 07 


19.067 45 
19.813 83 
19.838 42 


96.714 86 

96.750 70 

96.751 72 


1.192 415 
1.377 572 
1.132 323 


2.998 982 
4.088 539 
5.084 082 


5.325 360 
5.868 629 
6.735 687 


N(j max = 4) 


12.764 82 


20.040 74 


96.752 15 


1.579 273 


5.316 872 


6.317 204 



TABLE IV: The convergence of eigenenergies £t of Eq. (23) vs order A: max of approximation of 
effective potentials from (14) for j max = 4 and q max = 60 basis functions at 7^ = 0. Fast and slow 



variables xj = p and x s = z (prolate SQD and spherical QD), number of nodes i = (n pp ,n zp ). 



h 


c = 2.5, a = 0.5 


c = 2.5, a = 2.5 




(0,0) 


(0,2) 


(1,0) 


(0,0) 


(0,2) 


(1,0) 


8 

12 
20 


25.179 14 
25.199 62 
25.201 16 


34.076 77 
34.468 84 
34.522 02 


126.445 9 
126.456 
126.456 1 


1.471 911 
1.536 121 
1.563 492 


4.270 174 
4.716 984 
5.182 198 


5.614 892 
6.188 144 
6.266 533 


N(j max = 4) 


25.201 21 


34.525 12 


126.456 1 


1.579 239 


5.316 732 


6.317 058 



shows that for PSQD we have upper estimate and monotonic convergence with increasing 
&max to the numerical results at j max = 4. Similar behavior is observed for OSQD, however 
the accuracy of approximation of the effective potentials is worse, especially for the lowest 
effective potential %' — 1, corresponding to the ground state of the fast subsystem, because 
the upper estimates are violated. These Tables show also that such expansions have faster 
convergence for strongly oblate or prolate spheroidal QDs than for spherical ones. 
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TABLE V: The convergence of eigenenergies E% of Eq. (23) vs order /c max of approximation of 
effective potentials from (14) for j max = 4 and q max = 60 basis functions at = —10. Fast and 



slow variables Xf = p and x s = z (prolate SQD and spherical QD), number of nodes i = (n pp ,n 



zpj- 



k 

"'max 


c= 2.5, a = 0.5, 7 F = -10 


c= 2.5, a = 2.5, j F = -10 




(0,0) 


(0,2) 


(1,0) 


(0,0) 


(0,2) 


(1,0) 


8 

12 
20 


20.221 65 
20.607 33 
20.658 46 


30.913 36 
32.375 40 
32.674 45 


125.306 2 

125.331 6 

125.332 2 


-19.673 98 
-15.348 50 
-12.194 45 


-5.378 707 
-6.881 266 
-2.204 160 


-1.784 110 
-2.605 091 
-1.336 853 


N(j max = 4) 


20.66203 


32.708 77 


125.332 2 


-10.844 02 


-1.511 063 


1.129 039 



IV. PTRS FOR BVP FOR OSQD IN ELECTRIC FIELD BY FAST VARIABLES 



To have an analytic representation of the matrix elements (11) for small tjt, one can use 
Vf(xf,x s ) = 2^fz, Vf s (xf,x s ) = as potentials for OSQD instead of the potentials (12) 
introduced in Section 2.1. Then we arrive at the Sturm-Lioville problem for the OSQD in 
fast variable expressed in the form 

d 2 



dz 2 

(B^Bjifi)) 



ez-E J (p)^B 3 (z;p)=0, 



(24) 



L(p)/2 



L(p)/2 



Bi(z;p)Bj(z;p)dz = 5, 



in 



where e = 7^ is the electric field strength considered here as a formal parameter of the PT, 
implying a small interval p e (0, L(p) = 2ca/1 — p 2 /a 2 ) of the scalar product (Bi(p)\Bj(p)) . 
The solutions Bf\z; p) and Ef\p) of the unperturbed equation (at e = 0) have the form 



{B°\z;p),E°\p)} 



{B s j (z;p),E s j (p)}, for even j = 2, 4, 
{Bj(z;p),Ej(p)}, for odd j = 1,3,. 



(25) 



where 



Bfap) = y/2/L(p) sm(Trjz/L(p)), Bfap) = y/2/L(p) caa(njz/L(j>)), 
E°(p) = (nj/L(p)) 2 , Ej{p) = (nj/L(p)) 2 . 

We seek for the eigenfunctions BAz\ p) and the eigenvalues EAp) in the form of power 



expansions 



B 3 (z;p) = XVsf^p), E s (p) = J^e k Ef(p). 



(26) 



k=0 



k=0 
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Substituting Eq. (26) into Eqs. (24) and equating the coefficients at the same powers of e, 



we arrive at the system of inhomogeneous differential equations with respect to corrections 

3? 

d 2 



E {k) and Bf\z;p): 



dz 2 



Ef\p) Bf\z-p) = (z + Ef\p)) Bf- l \z-p) + Y J E { ?\p)Bf-v\z-p) 1 (27) 



p=2 



52(Bto(p)\B< k -»(p)) = 0. 

p=0 

In each k-th order of the perturbation theory (PT) the solutions becoming zero at the 
boundary points (z = ±L(p)/2 are sought in the form 

^ma x / j \ "max 2 , . 

ZB°(z;p)Sl k) z» + (z 2 -(L(p)/2) 2 ) E Bj(z; p)Cl%z", j = 2, 4, 



E B j( z 'i p)Cf ) z v + (z 2 - (L(p)/2) 2 ) E B*(z;p)Sl%z», j = l,3,... 

!/=0 Z^ = 



Substituting Eq. (28) into the corresponding equation (27) of the fc-th order of the PT, and 
extracting the coefficients at Bj(z; p)z v and Bj(z; p)z v , v = 0, i/ max , we arrive at the set 
of algebraic equations with respect to unknowns E^ k \p), Su and £ , for even j: 



-(-1)> + l)(L(p)/2)njCl% - (,/ + 2)(i/ + 1)4? 2 + (-1) J 2(^ + 1)^^ 
-^{p)^ ~ SO? -Y,Ef\p)Sr p) - Ef\p)5 vfi = 0, 

p=2 

+(1/ + + 2)(L(p)/2) 2 ^ fc + ) 4 - («/ + 2)0/ + l)Cj? a - (-1)^2(1/ + 1)^! 
-^(p)^ - {L{p)/2fC% X) ) - d£? + (L(p)/2fd k + ^) 

-Y,Ef\p){C^ - (L(p)/2M k + 7 ] ) = 0. 



p=2 



For odd j the same unknowns are calculated using the equations (29) (29) with the replace- 
ment S^. The unknowns ^ for even j and Sq for odd j are determined from the 
respective conditions: 

k 



EE^'^^wi^'i^w) ( 29 ) 

p=0 

+[(C7W5^ + + {L{p)l2)\d$ x &$ + 4 p i^;r ) )](^(p)i^ + l^( P )) 

+[(C«<#~ P) ) - 2(L(p)/2) 2 C« 1 C<#> + (L(p)/2)*C$ 3 C$$ }(B°(p)\z^'\B°(p))) , 
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and Sq for odd j is calculated from the equation (29) with the replacement 
This algorithm was implemented using the Maple environment. The run was performed 
until the maximal order of the PT k max = 8. Below we present the first few coefficients of 
the eigenvalue expansion, truncated by the terms proportional to e 6 = 7^ 

p * 2 j 2 , mP))\^3 2 ~ 15) , , (L(p)) 10 (1980 - 21Q 7 rY + 7 ry) 
^ (p) = WW 1 + + 2304^0,10 e ' ( 3 °) 



the eigenfunctions truncated by the terms proportional to e 2 = 7^ 



B-^P) + (-ffl^M + ^-yW~> ) £ , j = 2, 4, ... 



, . , (L(p))a*B?(*;p) L(p)( 2 2 _(£(p)/2) 2 )B|(*;p)\ . 

+ t- 2 e, J = 1,3,... 



3 v ' r 7 V 4-7r 2 j 2 47rj 

and the diagonal effective potentials, truncated by the terms proportional to e 6 = 7^ 



V dp y V 12 (^(p)) 2 576vr 6 / 



L(p) 10 (3510000 - 389880vr 2 j 2 + 3321vr 4 j 4 + 13vr 6 j 6 ) 4 . 
1 oTdo„i2„-i2 e • 



fdL(p) V / tt 2 j 2 + 3 (L(p)) 4 (-2880 + 2587r 2 j 2 + 7tt 4 j 4 ) 2 
I ' ■ ' ^ v 12(L(p))2 4 
89880tt 2 j 2 
276487T 12 j 12 



V. THE PTRS IN THE DIAGONAL ADIABATIC APPROXIMATION 

The desired solutions ofthe original 2D BVP ^ are determined by the diagonal approx- 
imation of the Kantorovich expansion ([7]) at fixed m 

The diagonal approximation of the BVP ([9]) and (10) in the slow variable has the form 

]_d d d + m _ \ = Q 

rr>LL si ry* fl K* r Y* I 

JU g \JjJU g \JjJU cf JU g J 

and the eigenfunctions satisfy the orthonormalization conditions on the semiaxis [x™ in = 
0, x™ ax = 00) at d — 1 for the OSQD and on the axis (x™ in = -00, x™ ax = 00) at d = for 
the OSQD 

„max 
/ x s 

XrA X s))Xi,n'{Xs){Xs) d dx s = 8 nn < . (33) 



Here Vi(x s ) = V s {x s ) + Ei(x s ) + DHu(x s ), where the parameter D is D = for the crude 
adiabatic approximation and D = 1 for the adiabatic approximation; V s (x s ) = 0, Ei(x s ) 



and Ha(x s ), Eqs. (30)-(31), for OSQD and K(x s ) = 27^,2, Ei(x s ) and Ha(x s ), Eq. (13) 
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FIG. 3: Three potential functions Vi(x s ) for oblate x s = p and prolate x s = z spheroids and their 
power expansions till sixth order with account of adiabatic frequencies Ui and lower bound shifts 
V {0) 

for PSQD; £i- n are the eigenenergies of a lower part of thespectrum £ i;0 < < ... < £ i;n 
enumerated in the ascending order by the number of nodes n = 0, 1, 2, ... of the eigenfunctions 
Xi-,n{x.s) at fixed adiabatic quantum numbers i = n Q for OSQD and i = n p for PSQD. The 
potential function Vi(x s ) is expanded in powers of the small parameter e 



V?-\x s ) = + + n~ 2 Yl;Z V i j) ^ 3 - ( 34 ) 

For OSQD at the values of the parameters d = 1, e = c~ 2 , k = 2, m = m the coefficients 



are determined by Taylor expansion of the effective potentials (|30h , (|3 1 h in the vicinity 



of the equilibrium point x s = 0. With the accuracy up to order of 0(7^) the coefficients 
V} j) and ujf are expressed as: 



2^2 



c 4 (vr 2 n 2 - 15) 4 4c 10 (7r 4 n 4 o - 2107r 2 n 2 + 1980) 



^ " 4c 2 +7f 3vr 4 n 4 +7f Ovrionio ' (d5j 



o 



7r 2 n 2 3 + 7r 2 n 2 2 / 8c 4 (7r 2 n 2 - 15) 4c 6 (77r 4 n 4 + 258vr 2 n 2 - 2880) 



i ~ (ac) 2 a 4 ~ lF \ 3aH 4 n 4 OuWn" 

+tI (- • ",.r;;r:,; " " ••' ■ 



( 80c 10 (% 4 n 4 - 210vr 2 n 2 + 1980) 16c 12 (137r 6 n6 + 332l% 4 n 4 - 389880vr 2 n 2 - 3510000) \ 
V 9A%i + 27a% 12 ni 2 ) 



T/ (i) _ ^ . n 3 + 7 r 2 ng 2 / 4c 4 (7r 2 n 2 -15) ^ n 4c 6 (77r 4 n 4 + 258^ - 2880) 

^ ~ VM 2 « 4 7F V 3a%%4 d < 2 - ^ g^v^e d 

4 /I6c 10 (vr 4 n 4 - 210^ + 1980)^ nx 1n fe lK 6 i4 <J«. 



;2 



+1f[ ° Q 4 10 10 ° ^(105 J2 -10^ + 5 , 

\ 9a 4 7r lu n i , u a 1 a 4 a b 

16c 12 (137r 6 ^ + 332l7r 4 n 4 - 38988Q 7 r 2 n 2 , - 3510000) L kn -, ;+2 

+ 27a6vr 12 ni 2 1 ° l2 + a 2 a 4 + a 6 ' 1 1 * 
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For PSQD at the values of the parameters d = 0, e = 1, k — 1, m = the coefficients are 



sought in the form of a Taylor expansion in powers of x s = (x s — Xq) and 7^ of the effective 
potentials Vi{x S) ^p) = Ej(x s ) + DHjj(x s ) + Jf^s, Eq. (13). The expansion coefficients 
x = ^ k T 2 k+il F k+l are sought from the equilibrium condition dVl ^' lF ' ) = at fixed 

s x s =xo 

7jt. With the accuracy up to 0{^\) the coefficients and uf are expressed as: 

= -2r lT | - 2r 37 ^ + < iH /(a 2 ) + D(l + < iH )/(3c 4 ) 
+7^ 1 2 (< j|m| /(a 2 c 2 ) + D(l + a 2 p>|m| )2/(3c 4 )) 

+7Fn« i | m ,(r 1 3 + 2c 2 r 3 )/(a 2 c 4 ) + 2D(1 + a 2 npM )(r? + r 3 c 2 )/(3c 6 ), 



CO; 



a 



n p ,\m\ 



/{a 2 c 



)/(3c 4 



(36) 



+ 7 |r 2 6(a 2 pi | m |/(a 2 c 4 ) + D(l + a 2 pj|m| )2/(3c 6 )) 

+7ki«,| m |(15r 1 3 + 12c 2 r 3 )/(a 2 c 6 ) + D(l + < fH )(157f + 8r 3 c 2 )/(c 8 ) 
= (+27FTx((i + l)<, H /(a 2 c 2 * +2 ) + D(i + 1) 2 (1 + <, H )/(3c 2i + 4 )) 

+2 T 3 (i + l)(a 2 pjH ((2z 2 + 7z + 6)r 3 + 3r 3 c 2 )/(3a 2 c 2i+4 ) 
2 

a p ,|m| 



+£>(! + <. m .)((2i s 7f + ll* 2 + 20* + 12)7f + 3(i + l)r 3 c 2 )/(9c 2 *+ 6 ) ) +1 



+ a. 



n p ,|m| 
2 ^2 



/(aV l+2 ) + L>(l + < |H )(z + l)/(3c A+4 ) 



+ 7 2 r 1 2 (, + 2)(2< + 3)(< H (aV» +4 ) + £>(1 + < )|m| )(i + 2)/(3c 2 ^)) 
+7^n(i + 2)(2i + 3)(a 2 p;|m| ((2* 2 + llz + 15)r x 3 + 12c 2 r 3 )/(6a 2 c 2l+6 ) 

+D{\ + a 2 pj|m| )((2z 3 + 17z 2 + 48i + 45)r 3 + 12(i + 2)r 3 c 2 )/(18c 2i+8 

where T2/c+i is determined from the condition that the coefficient at x s is zero: 
3a 2 c 4 \<i"< -'■■(/)< ... .1 n 



2i+2 



3c 2 a? 1 ™ 1 + J Da 2 (l + <, H ; 



54aV°(3c 2 « 2 piH +2 J Da 2 (l 



n p ,\m\ 



)) 



t I 

n p ,|m| 



(3c 2 a 2 , , + L>a 2 (l + a 2 , , 

V n p ,\m\ V n p ,|m| 



)) 4 



In Fig [3] we show three potential functions V^(x a ) for oblate x s = p and prolate x s = z 
spheroids and the convergence of the corresponding power expansions till the sixth order 
with account of adiabatic frequencies ajj and lower bound shifts 



We choose the unperturbed operators of Eq. (32) at e = in the expansion (34) in 



the form ( 16 )— ( 19 ) with the eigenvalues and the basis functions of 2D- and ID- oscillators 
given in Section 3 with respect to the scaled coordinate x, x s = \j2x/ui and x s = x/y/uo^ 



where the adiabatic frequencies Ui are defined by Eqs. (35) and p6J) (at fixed i' = n + 1), 



respectively. According to (34), we seek for the eigenfunctions Xi;n( x s) an d the eigenvalues 
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(k) (k) 

Si- n in the form of expansions in powers of e with unknowns and E n , omitting the 
notation m for brevity: 

^max 

X,n(x s ) = < 0) + ^$i fc) (x s )e fc 

k=l 

^max I ^max 

= + E4!? = + ^ 0) + E E ^ h 



(37) 
(38) 



fc=0 



fc=i 



Substituting the expansions (34), (37) and (38) into Eq. (32) and equating the terms with the 



same power of the parameter e, we arrive at the recurrence set of inhomogeneous equations 
of the PT with respect to the unknowns E n k ^ and &n\x): 



L{n)*S\x) = = /(°)(x), 



(39) 



fc-i 

L{n)$W(x) = J2( E t P) ~ Vt p) )^\x) ee /«(*), k > 1. 

p=0 



with the initial conditions (16) and (18) for OSQD and PSQD, respectively. The solution 



of this problem is implemented in four steps. 



Applying the relations (17) and (19), we expand the right-hand side f( k '(x) and the 
solutions $( fe )(x) of Eqs. (39) over the basis of normalized states $^ s (x), Eqs. (16) and 



(18): 



(40) 



Then a recurrent set of linear algebraic equations for unknown coefficients bi k ^ and corrections 
E^ is obtained 



> h (k) _ ,(fc) Q _ 
J J s w 5 ° ''man • • • 5 °maxi 



(41) 



where s' — s for OSQD and s' = 2s for PSQD. These equations are solved sequentially for 



h — 1 9 

ft j. , , . . . , ru max . 



;( fc ) = ^^( fc ); 6( fc ) = /W/ s ', s = - s 



The initial conditions for this procedure are 



'1.x; • • • i Smaxi S ^ 0. 



(42) 



6< 0) = 5 s0 , E<® = (n + (\m\ + l)/2) or = (n + l)/2). 
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To obtain the normalized wave function &j(x) up to the k-th order, the coefficients b 
are determined by the following relation: 

1 Sma. 



(k) 



■,(*) 



2(0|0) 



(43) 



p=l s'- 



The above scheme implemented in Maple was applied to the evaluations of solutions in 
the analytical form up to the order k max = 6 of the PTRS. The first four nonzero coefficients 



for the energy (38) in the analytic form, truncated by the terms proportional to the sixth 
power of the electric field strength, 7J1, in the crude adiabatic approximation (CAA) take 
the form: 

1) For OSQD in terms of minor c and major a semiaxes; the set of adiabatic quantum 
numbers [m, n Q = n zo + 1, n po ] 



K (0) 

ItlQ 

£ (0) 



£(1) 



ir 2 n 2 



Ac 2 



- + r F 



c 4 {jt 2 n 2 - 15) 4 4c 10 (7r 4 n 4 o - 210vr 2 n 2 + 1980) 



3ir 4 n 4 



97r 10 ni 



(44) 



TTUo 

ac 
1 



2 Ac 5 (ir 2 n 2 - 15) 4 8c 11 (2Tr 4 n 4 - 360vr 2 n 2 + 3375) 
7f ir~ t-r If 



Or 



+ 1F 



3air 5 nl 
Ac % (ix 2 n 2 - 15) 



3air 11 nl 1 



(2n po + \m\ + 1), 



4 16c 12 (77r 4 n 4 - IIIOttX + 10350) 

+ 1f 



3ir 12 n} 2 a 2 



x(2 + 6n po + 3\m\ + 6n 2 po + \m\ 2 + Qn po \m\) 

3c 



4% 00 = ( 2n P o + H + i 



2na 3 n f 



(2 + 2n po + \m\ + 2n 2 + 2n po \m\) 



-If 



-If 



2 2C 7 (7T 2 ^ 



15) 



(54 + 118n po + 16|m| 2 + 59|m| + HSn 2 po + 118n po |m 



3n 7 n 7 a 3 

4c 13 (1874vr 4 n 4 . - 273120vr 2 n 2 + 2536425) 



97r 13 nl 3 a 3 



(2n po + \m\ + 2n 2 + 2n po \m\) 



224c 13 (8vr 4 ^ - lMOvrX + 10575) 



m + 



8c 13 (326tt 4 ™ 4 - 48480vr 2 n 2 + 450675) \ 



3n 13 nl 3 a 3 

2) For PSQD in terms of minor a and major c semiaxes, the set of adiabatic quantum 
numbers [m,n p = n pp + 1, n zp ) and positive zeros a np ,\m\ of the Bessel functions of the first 
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kind 



V (o) 



£(P) 

Tlpyf 



a 



n p ,\m\ 
a n p ,\m\ 

ac 



If 



a 2 c 2 

A 2 

n p ,|m| 



+ 7^ 



a 6 c 4 



16 <,h' 



(45) 



3a 3 c 
4a 3 , , 

n p ,\m\ 



-If 



9a' c 



7^3 



16a 7 , i 

n p ,\m\ 



(2n zp + l), 



C(2) 



3 2 27a 4 
4^ + 7F 16a 4 
3a 



If 



105a 8 c 2 
64a 1 



(2n 2 + 2n zp + 1), 



"■" 16c 3 a ■ ,( 2n ^ + 1 )« + n ^ + 3 ) 



+7l 



-7^ 



*n p ,|m| 

5a 5 



64c <,M 
45a 9 c 



(2n 2p + l)(25n 2 p + 25n zp + 51) 



(2n zp + l)(23n^ + 23n zp + 37) 



4a n P ,\m\ 

3a 8 c 2 



(30^ + 30^ + 11 



5a 



(30n 2 +30n zp + 11) 



n p ,|m| 



In Tables |VI| and |VII| we demonstrate how the approximate eigenvalues in the lower part 
of spectrum for OSQD and PSQD at m = and 7^ = converge to the values calculated 
numerically with required accuracy in the crude adiabatic approximation with increasing of 
the PT order k. The accuracy was from 8 to 5 digits at n zo = 0, from 10 to 8 digits at 
n zo = 2, from 6 to 4 digits at n pp = 0, and from 8 to 7 digits at n pp = 1, respectively. Note, 
that the difference between the adiabatic shift and the eigenvalues £^ n = + Sj® in 
the zero order k = of the PT is small, but increases with growing n po and n zp for OSQD 
and PSQD, respectively. The shifts give the main contribution and provide the lower 
adiabatic estimate of each set of eigenvalues, generated by the perturbed harmonic oscillator 
terms with adiabatic frequency Wj. From Tables VI and VII] one can see that with increasing 
quantum numbers n zo (or n pp ), related to the fast variable, the accuracy of approximation 
of the lower part of the spectrum is increasing. This is because the accuracy of the Taylor 



approximations of potential function (|34|) in Eq. (|32|) is improved with increasing the number 

i = n zo + 1 > 2 (or i 



n pp + 1 > 2), which is demonstrated in Fig. |3j 
In Figs. |4]and|5]we show the eigenvalues S of the lower part of the spectrum of oblate and 
prolate QDs versus the electric field strength within small (left panels) and large (right pan- 
els) intervals of jp, calculated in the crude adiabatic approximation (solid and dashed lines) 
to compare them with the numerical results (dotted lines). One can see that the eigenvalues 
calculated using the PT (solid and dashed lines), corresponding to the eigenfunctions with 
smaller number of nodes along the electric field (i.e., with smaller n zo for OSQD and n zp 



24 



TABLE VI: Convergence of eigenvalues = + 42,n po for o&tote spheroid c = 0.5, 

a = 5 vs PT order k max at 7^ = . First line * notes adiabatic shift Vn®) ,n po - Last line are results 
of numerical calculations (Num). 



I 

"-max 


Tlzo — 0,71p O — 


n zo = 0,n po = 


1 n zo = 0,n po = 


2 n zo — 0,fip O — 


3 n zo — 0,rip O — 


4 




11.12624146 


13.63951558 


16.15278970 


18.66606383 


21.17933795 







11.20624146 


14.19951558 


17.67278970 


21.oz60o383 


26.05933795 




1 


11.21006118 


i4t.zoooyouo 


1 f .OUD4 f yoo 


91 Q7'?R I: ;S99 

zi.y / oDOozz 


9fi 781 9R/I77 
ZD. 1 O1Z04 1 1 




2 


11.21026382 




1 7 8095485Q 


91 951 00981 


96 710Q4787 




3 


11.21028027 


1 4 9^441 79'} 


1 7 80949765 


91 94Q08610 


96 7095691 5 




4 


11.21028227 


1 4 9^44' : i790 


1 7 809651 Q5 


91 95065951 


96 7095Q1 63 




5 


11.21028259 


1 4 9^444049 


1 7 80964785 


91 950599Q1 


96 70875037 




Num 11.21028268 


1 4 9*34441 47 


1 7 809^0^ 




9fi 708^7797 

ZU. 1 UOJ / / Z 1 




h 

"rnax 




TL zo — 2,7ip — 






3 ?1-2o — 2,71p — 


4 


* 


92.59635079 


100.1361731 


107.6759955 


115.2158178 


122.7556402 







92.67635079 


100.6961731 


109.1959955 


118.1758178 


127.6356402 




1 


92.67762403 


100.7076323 


109.2405589 


118.2916826 


127.8762825 




2 


92.67764654 


100.7076818 


109.2401221 


118.2891654 


127.8684695 




3 


92.67764715 


100.7076847 


109.2401176 


118.2890944 


127.8681589 




4 


92.67764718 


100.7076850 


109.2401203 


118.2891137 


127.8682457 




5 


92.67764718 


100.7076850 


109.2401203 


118.2891132 


127.8682422 




Num 92.67764718 


100.7076850 


109.2401204 


118.2891132 


127.8682419 





for PSQD) and with greater number of nodes across the electric field (i.e., with greater n po 
for OSQD and n pp for PSQD), provide better approximation of the eigenvalues, calculated 
numerically with required accuracy (dotted lines). This property follows from the fact that 
such functions have better localization in the vicinity of the plane, passing through the QD 
center transverse to the electric field, i.e., in the region with minimal contribution of the 
electric field potential to the Hamiltonian of the system. As shown in the right panels of 
Figs. |4]and|5j the differences between the egienvalues, calculated using the PT and the nu- 
merical method, increase faster in a smaller interval of 7^ for larger PSQD than for smaller 
OSQD, the size being measured along the direction of the electric field. 
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TABLE VII: Convergence of eigenvalues 4^*1 = V Zl + ^Ifeo* £ n%,n zp for prolate spheroid 
c = 2.5, a = 0.5 vs PT order A; max at 7^ = . First line * notes adiabatic shift Vn p ,n zp - Last line 
are results of numerical calculations (Num) . 



^max 


Tipp — 0,n Z p — 


Tlpp — 0,7i.2:p — 


2 Tlpp — 0iTl Z p — 


4 Tipp — : n Z p — 


6 ripp = 0,n zp = 8 


* 


25.05660430 


32.75204608 


40.44748787 


48.14292965 


55.83837144 





25.17660430 


34.31204608 


45.36748787 


58.34292965 


73.23837144 


1 


25.18408925 


34.42432034 


45.88394944 


59.80249498 


76.41947535 


2 


25.18465987 


34.42810718 


45.87103269 


59.69485535 


76.04779441 


3 


25.18472054 


34.42867746 


45.87114189 


59.68460238 


75.99436976 


4 


25.18472960 


34.42880826 


45.87257549 


59.69618640 


76.05191800 


5 


25.18473139 


34.42883580 


45.87259458 


59.69511288 


76.04351256 


Num 25.18472985 


34.42884694 


45.87265876 


59.69512314 


76.04210082 


^max 


Tipp — \,n Z p — 


Tlpp — Ij^zp — 


2 Tipp — \,n Z p — 


4 Tipp — l,n Z p — 


6 ripp = l,n zp = 8 


* 


126.3011119 


143.9653618 


161.6296118 


179.2938617 


196.9581117 





126.4211119 


145.5253618 


166.5496118 


189.4938617 


214.3581117 


1 


126.4243727 


145.5742742 


166.7746086 


19 0.1297223 


215.7439616 


2 


126.4244810 


145.5749929 


166.7721571 


190.1092932 


215.6734198 


3 


126.4244860 


145.5750400 


166.7721661 


190.1084455 


215.6690025 


4 


126.4244863 


145.5750447 


166.7722178 


190.1088627 


215.6710754 


5 


126.4244864 


145.5750452 


166.7722181 


190.1088459 


215.6709435 


Num 126.4244896 


145.5750487 


166.7722220 


190.1088484 


215.6709278 



The range of the parameter values, for which the PT algorithms are valid, was estimated 
by means of numerical calculations using the KANTBP program [30] , as well as the condition 
that the mean value of the slow variable is smaller than the size of the major axis of OSQD 
or PSQD, i.e., p < a or z < c, or known estimates of the distribution of nodes of Laguerre 
or Hermite polynomials [37]. To calculate also the approximate eigenfunctions of the lower 
part of the spectrum n = 0,...,n max with required numbers n of nodes in the interval 
p E (0, a) (or z G (— c, c)) for OSQD (or PSQD), one should choose such value of parameter 
a = 2x 1 Ui (or c = x /y/cJi), that outside this interval x G (x = An + 2\m\ + 2, 00) (or 
\x\ G (x = (2n + l) 1/2 ,oo)) the Laguerre (or Hermite) polynomials have no nodes. As an 
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FIG. 4: Dependence of eigenenergies £ (in units of E e ) of lower part of spectrum of electronic 
states of OSQDs (a = 2.5, c = 0.5) at m = on electric field strength 7^ (in units of Fq ): Solid and 
dashed lines are eigenenergies calculated by PTRS till 5 order in crude adiabatic approximation: 
seven solid lines (n zo = 0, n po = 0,1,..., 6) and four dashed lines (n zo = 1, n po = 0,1,2,3) are 
shown on left panel in interval 7^ G (0, 1); the first three of solid lines (n zo = 1, n po = 0,1,2) and 
the first two of dashed lines (n zo = 1 , n po = 0, 1) are shown on lower-right and upper-right panel, 
respectively, in bigger intervals jp G (0,20) and jf £ (0, 10). Numerical solutions of Eqs. Q at 
jmax = 4 are shown by points. 

example, in Fig. [2] we show contour plots in zx and xz plane of the first four eigenf unctions 
of OSQD and PSQD, respectively, that have a required number of nodes (crossings of the 
function plot with zero plane) in the interval p e (0, a) and z G (— c, c) at the values c = 0.5, 
a = 5, c = 2.5, a = 0.5. One can see that the asymmetry with respect to z-axis of the 
eigenf unctions of PSQD is greater than that of OSQD, because the variation of well depth 
of PSQD is greater than of OSQD. 
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FIG. 5: Dependence of eigenenergies £ (in units of E e ) of lower part of spectrum of electronic 
states of PSQDs (a = 0.5, c = 2.5) at m = on electric field strength 7^ (in units of Fq): Solid and 
dashed lines are eigenenergies calculated by PTRS till 5 order in crude adiabatic approximation: 
fifteen solid lines (n pp = 0, n zp = 0, 1, 14) and three dashed lines (n pp = 1 , n zp = 0, 1, 2) are 
shown on left panel in interval G (0, 0.1); the first six of solid lines (n pp = 0, n zp = 0, 1, 5) are 
shown on lower-right panel and the dashed lines (n pp = 1 , n zp = 0, 1, 2) are shown on upper-right 
panel in bigger interval £ (0,3). Numerical solutions of Eqs. Q at j max = 4 are shown by 
points. 

VI. ABSORPTION COEFFICIENT FOR AN ENSEMBLE OF QDS 

One can use the differences in the energy spectra to verify the considered models of QDs by 
calculating the absorption coefficient K(u ph , a, c, ) of an ensemble of identical semiconductor 
QDs [14J. Since we do not discuss exciton effects in the present paper, the absorption 
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coefficient may be approximately expressed as 

K(oo p \ a, 5,u) = J2 Kvy(£o p \ a,d,u) = AJ^ I»y(u)5(Hoj ph - W vv >), (46) 



vy vy 

2 

%(r; a, c, F, /x e )^/(f ; a, c, F, fi h )d 3 r 



where A is proportional to the square of the matrix element in the Bloch decomposition, 
^l(u) and \&^/ are the eigenfunctions of the electron (e) and the heavy hole (h), Ef, and 
E„, are the energy eigenvalues for the electron (e) and the heavy hole (h), depending on the 
semiaxis size c, a for OSQD (or a, c for PSQD) and the adiabatic set of quantum numbers 
v = [n zo ,n po ,m\ and v' = [n' zo ,n po >,m'\ {y = [n pp ,n zp ,m] and v' = [n' pp , n' zp , m']), where 
m! = —m, Eg is the band gap width in the bulk semiconductor, u ph is the incident light 
frequency, W vv > = E g + E^(a,c) + E^,(a,c) is the inter-band transition energy for which 



K(u ph ) has the maximal value. We rewrite the expression (46) in the terms of frequency 
shift of the incident light Au ph / (2tt) = (hu ph — E g ) / (2irh) corresponding to the inter-band 
transition energy shift AW vv i = W vv > — E g = E%(a, c) + E%,(a, c) for which K(A£j ph ) has the 
maximal value, using dimensionless variables in the reduced atomic units 

K(Au p ,a,c) = AE g ^ I u y(u)5[f u y(u)\, f v y[u) = X t — | (47) 

i>y & 

where the parameter u will be defined below, Ai = (hu ph — E g )/E g is the energy of the optical 
interband transitions scaled to E g , 2E g = E g /E R is the dimensionless band gap width. 

For GaAs the functions f^ e (u) describing the (h — > e) interband transitions have the 
form 

fuy e (u) = Xi- (2E g )- 1 (2Et(a,c, lF ) + 2E e v ,(a,c,-(fi h /fi e ) lF )(fi e / f i h )), (48) 

where \i e = 0.067m and = Hhh = 0.558m are the masses of electron and hole, respec- 
tively, Eg = 1430 meV is the band gap width and k = 13.18 is the dc permittivity and 
Ex = e 2 /(2na e B ) = 5.275 meV, a e B = h 2 K/(fi e e 2 ) = 104 A, E h R = e 2 /(2Ka h B ) = 49 meV, 
a% = h 2 K/(fi h e 2 ) = 15A, 2j F = F/F*, F* = E e R /{ea B ) = e/(2/t(a^) 2 ) = 5.04kV/cm. 

For InSb the dispersion law for heavy holes (hh) is parabolic while for electrons (e) and 
light holes (Ih) it is non-parabolic and may be described by the Kane model [HI [19l [22] at 



29 



7f = 0. The energy values in our notation are: 

2E™(InSb) = 2E*,(a, c), (49) 



2E e u {InSb) = 2E l v h (InSb) = -E g /2 + jE 2 g /A + E g (2E^(a, c)). (50) 



As follows from Eqs. (49) and (50), to determine the energy spectrum and the wave 
function of the light hole and the electron one should solve the Klein-Gordon equation ([39, 
140]). while for heavy hole the Schrodinger equation is applicable. The functions f^^ e {u) 
and f„^ e (u) describing the (hh — > e) and the (Ih — > e) interband transitions have the forms 

= Ai - (1/2 + y/l/A + (2EZ(a, c)/{2E g )) + (2^)- 1 2^,(a, c)(^//x A )), (51) 

ft^» = Ai - 2(y/l/4 + (2EZ(a,c)/(2E g ))), (52) 

where /x e = = 0.15m and /i^ = = 0.5m are the masses of electron, light and heavy 
holes, respectively, E g = 180 meV is the band gap width, k = 16 is the dc permittivity, 
and E e R = E% = e 2 /(2na e B ) = 7.972 meV, a% = a% = h 2 n/(fi e e 2 ) = 56.44A, E\ = Ef = 
e 2 /(2Ka B h ) = 26.57 meV, a h B = af = h 2 n/(fx h e 2 ) = 16.93A. 

For both electron and hole carriers the dimensionless energies 2EI = Ef,/E R and 
2E^(/j lfl / fMe) = E^/E R are expressed in the same reduced atomic units E R , and the over- 



lap integral (46) between the eigenfunctions, corresponding to E^p) and E^{^p) 



{fie/ ^h)E^(-(fi h / fie^p), takes the form 



h.v' (u) 



2 

( a s) 3 ^( r ; a > c, Jf, fi e )*£/(r; a, c, -(jih/He)lfF, He)d 3 r . 

Now consider an ensemble of OSQDs (or PSQDs), differing in the minor semiaxis values 
c = u Q c (or a = u p a), determined by the random parameter u = u Q (or u = u p ). The 
corresponding minor semiaxis mean value is c at fixed major semiaxis a (or a at fixed major 
semiaxis c), and the appropriate distribution function is P(u Q ) (or P(u p )). Commonly, in 
this case the normalized Lifshits-Slezov distribution function [15] is used: 

P{u) = {3 V exp(-l/(l - 2u/3))/2 5/3 /(u + 3) 7/3 /(3/2 - u) 11/3 , u G (0, 3/2); 0, otherwise} 

having conventional properties j P(u)du = 1, u — J uP{u)du = 1. The absorption coeffi- 
cients K°(u ph , a, c) or K p (u ph , a, c) of an ensemble of semiconductor OSQDs or PSQDs with 
different dimensions of minor semiaxes are expressed as 

K°(tu ph , a, c) = [ k(u ph ,a,c,u )P(u )du , K p (u ph ,d,d) = [ K(u ph ,d,d,u p )P(u p )du p . 



(53) 



(54) 
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Substituting (47) into (54) and taking into account the known properties of the <5-function, 
we arrive at the analytical expression for the absorption coefficient K(u ph ,a,c) of a system 
of semiconductor QDs with a distribution of random minor semiaxes: 



K(u ph ) 
K 



K v y{Q ph ) 
K 



K u y{Cj ph ) 



hy (u s 



df„y{u) 




du 


u=u s 



(55) 



where Kq = A l E g is the normalization factor, u s are the roots of the equation f u y(u s ) = 0. 

At 7jt = for IPBM we have the interband overlap I v y = &n po y p J>n zo y z J>m,-m' 
for OSQD, or I v y = (Ji + | m |(a npp+ ij m |)/Ji„| m |(a npp+ ij m |)) 2 5 nzpi „ ip 5 nppin ^5 m _ m / for PSQD, 
where a n +i t \ m \ is the positive root of the Bessel function, and the selection rules m = —m', 

n' zp [27], while at 7^ 7^ one should calculate 



n, 



n poi OT n pp 



n ppi 71 zp 



the interband overlap (53) in accordance with the selection rules m = —m', n po 



n P oi or 



n 



pv 



n 



ppi 



respectively. Note, that in the adiabatic limit and at small jf the contributions 



of non-diagonal matrix elements to the energy values are about 1% for IPBM of OSQD and 
PSQD; then in the Born-Oppenheimer approximation of the order b max for the AC we get 



,1 (u) 



fmax 

j=0 



(56) 



The coefficients of the expansion (56) for parabolic dispersion law for small 7f 7^ were 



constructed using the expansions (44) and (45) and at 7^ = they are given in |27J. In 



general case for the calculation f u y(u) by formula (48), (51), or (52) we used the eigenvalues 
E^(a,c) and E^;,(a,c) calculated numerically with given accuracy. After that we evaluated 



the coefficients of expansion like (56) by the method of least squares and by the polynomial 
interpolation in the case of parabolic and non-parabolic dispersion laws, respectively. Be- 
cause of monotonic behavior of function f v y(u) vs u in the case under consideration, we 
have only one root u s of the equation f u y(u s ) = 0, which was used in formula (55). 

For the Lifshits-Slezov distribution Figs. [6] and [7] display the total absorption coefficients 
K(u ph ) / 'Kq and the partial absorption coefficients K uu (u ph ) / 'K , that form the correspond- 
ing partial sum (55) over a fixed set of quantum numbers z/, v' at m = —vnl = 0. As a 



result of averaging (54) a series of curves with finite with and height are observed instead 



of a series of ^-functions. One can see that the summation over the quantum numbers 



n. 



+ 1 



1, 2, 3, 4, 5 (or n p 



n 



pp 



+ 1 = 1,2,3) enumerating the nodes of the wave 



function with respect to the fast variable gives the corresponding principal maxima of the 
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FIG. 6: Absorption coefficient K/Kq, Eq. ( 55 ), consisting of a sum of the first partial contributions 



vs the energy A = Ai of the optical interband transitions for the Lifshits-Slezov distribution, using 
the functions f^ e (u) for GaAs (h — > e) without electric field: for ensemble of OSQDs c = 0.5, 
a = 2.5 (left panel) and for ensemble of PSQDs a = 0.5, c = 2.5 (right panel). 





FIG. 7: The same as in Fig. [6j but in the presence of electric field 2^p = F/Fq. For comparison, 
the corresponding absorption coefficient without electric field is given by dashed line. 

total AC for the ensemble of QDs with distributed dimensions of minor semiaxis, while the 
summation over the quantum number n po = 0, 1, 2, 3, 8 (or n zp = 0, 1, 2, 15) that labels 
the nodes of the wave function with respect to the slow variable leads to the increase of 
amplitudes of these maxima and to secondary maxima arising in the case of sparer energy 
levels of IPBM of OSQDs (or PSQDs). 

In the regime of strong dimensional quantization the frequencies of the interband tran- 



sitions (h — > e) in GaAS between the levels n = l,n 



po 



0, m 



for OSQD or 



n.,. 



^,n zp = 0,m = for PSQD at the fixed values a = 2.5a e and c = 0.5a e for 
OSQD or a = 0.5a e and c = 2.5a e for PSQD, are equal to Aa)f 00 /(27r) = 16.9THz at 
7 F = and Au p 1 q /{2ti) = 15.9THz at 7 F = 10, or Au^ /(2tt) = 33.3THz at 7f = and 
Ao)^oo/(27r) = 31.5THz at jp = 2, where Aa}^Q /(27r) = (27r/i) _1 (H / ioo,ioo — E g ) corresponds 
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FIG. 8: The same as in Fig. rol but for InSb (hh — > e) interband transition. 




2.0 
1.5 
1 1.0 
0.5 










/"v 




f 

/ 








/ 




- '-:v~ ~.-vj?^3 




fff, 

/ 1 I 

I ' ' ' 1 ■ 





10 xi 



10 



A1 



15 



FIG. 9: The same as in Fig. |6j but for InSb (Ih — > e) interband transition. 

to the IR spectral region E], taking the band gap value {^Inh^Eg = 346 THz into ac- 
count. In Fig. [7] one can see the quantum-confined Stark effect that consist in the reduction 
of the absorption energy (light frequency) at the expense of lowering the energy of both (e) 
and (h) bound states due to the electric field effect. The total ACs at F ^ 0, shown by 
solid lines, qualitatively correspond to the total AC at F = 0, shown by dashed lines, but 
have lower magnitudes and smooth behavior, in spite of the additional contribution to the 



partial ACs of the overlap integral (53) from the interband transition n zo ^ n' zo or n zp ^ n' zp 
in OSQD or PSQD, also shown in Fig. [7j 

At the same parameters of the QDs the frequencies of the interband transitions (Ih — > e) 
in InSb are equal to Aa)f£ /(27r) = 68.5THz for OSQD or Au p 1 q /(2tt) = 87.2THz for PSQD, 
while the frequencies of the interband transitions (hh — > e) in InSb are equal to Ao)^/ (2tt) = 
78.6THz for OSQD or Au P q /(2tt) = 102THz for PSQD. These values correspond to the 
infrared spectral region with longer wavelength, similar to |22j, with the band gap value 
(2'n;K)~ 1 E g = 44THz taken into account. One can see that the behavior of total ACs for 
parabolic dispersion law for IPBM of InSb, shown in Fig. |8l is similar to that for GaAs (Fig. 
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[6]), while the behavior of AC for nonparabolic dispersion law, shown in Fig. [9j is essentially 
different. In particular, for OSQDs it grows faster with increasing Ai, while for PSQDs it 
goes to a plateau before starting to grow. Indeed, with increasing quantum numbers n po or 
n zp that characterize the excitation of slow motion, the maxima of partial ACs decrease for 
parabolic dispersion law, while for the nonparabolic one the maxima of partial ACs increase. 

With decreasing semiaxis the threshold energy increases, because the "effective" band 
gap width increases, which is a consequence of the dimensional quantization enhancement. 
Therefore, the above frequency is greater for PSQD than for OSQD, because the SQ im- 
plemented in two direction of the plane (x,y) is effectively larger than that in the direction 
of the z axis solely at similar values of semiaxes. Higher-accuracy calculations reveal an 
essential difference in the frequency behavior of the AC for interband transitions in systems 
of semiconductor OSQDs or PSQDs having a distribution of minor semiaxes, which can be 
used to verify the above models. 



VII. CONCLUSION 



The 3-D BVP for spheroidal quantum dots with respect to fast and slow variables of 
cylindrical coordinates was reduced by Kantorovich or adiabatic method to BVP for set of 
second-order differential equations (ODE) with effective potentials given in the analytic form 
with respect to the slow variable, using the basis function of fast variables, that depended 
on the slow variable as a parameter. Separation of variables of 3D BVP in spheroidal coordi- 
nates provides exact classification of energy eigenvalues by means of nodes of eigenfunctions 
which transforms exactly to an adiabatic classification of eigensolutions of a diagonal ap- 
proximation of ODE at small parameter, i.e. ratio of minor and major semiaxes of oblate 
or prolate spheroid. The effective potential of a crude diagonal adiabatic approximation 
(CDAA) of the ODE has been approximated by power expansions by slow variable. Energy 
eigenvalues and eigenfunctions of the BVP for CDAA were sought in expansions over eigen- 
functions of 2D or ID oscillator with adiabatic frequencies and power of small parameter 
by the PT. Required coefficients of these expansion were calculated in analytical form as 
polynomials of the sets of adiabatic quantum numbers. 

To specify the region of the model parameters, in which the PT asymptotic series are 
valid, we we compared the PT results with those of numerical calculations carried out with 
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required accuracy. The PT eigensolutions were used in analytic evaluation of the photoab- 
sorption coefficient for ensembles of oblate and prolate spheroidal QDs with given random 
distribution of small semiaxes without and with small values of external electric fields. In 



general case for calculation f u y(u) by formula (48), (51), or (52) we used eigenvalues calcu- 



lated numerically with given accuracy and we evaluated the coefficients of expansion like (56 ) 
by the method of least squares and by the polynomial interpolation in the case of parabolic 
and nonparabolic dispersion laws, respectively. Note, in the case of numerical calculations of 
the photoabsorption coefficient the required derivatives of eigenenergies and eigenfunctions 
with respect to a parameter, e.g., the small semiaxis, can be calculated also with the help 
of the numerical algorithms [291 135] . 

The elaborated methods, symbolic-numerical algorithms (SNAs) and programs [23W35] 
can be applied for solving the BVPs of discrete and continuous spectra of the Schrodinger- 
type equations and the analysis of spectral and optical characteristics of QWs, QWr's and 
QD's in external fields, as well as the spectra of models of deformed nuclei [41] . 

This work was partially supported by the RFBR Grants No 10-02-00200 and 11-01-00523. 
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